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ABSTRACT 


An analytic theory is developed for the time dependent 
magnetic fields inside the Moon and the diamagnetic cavity 
when the interplanetary electromagnetic field fluctuation 
propagates parallel to the cavity axis. The Moon model has an 
electrical conductivity which is an arbitrary function of radius. 
The lunar cavity is modelled by a nonconducting cylinder extend- 
ing infinitely far downstream. ..For frequencies less than about 
50 Hz, the cavity is a cylindrical waveguide below cutoff. 

Thus, cavity field perturbations due to the Moon do not propa- 
gate down the cavity, but are instead attenuated with distance 
downstream from the Moon. Far from the Moon, the cavity electro- 
magnetic field is a cylindrical TE mode propagating downstream 
with the same frequency and wavelength as the interplanetary 
field. Thus the magnetic field in the far downstream cavity has 
a component parallel to the cavity axis which is 90° out of phase 
with the incident magnetic field. The far cavity field is the 
result of a surface wave on the cylindrical boundary induced by 
the interplanetary field and moving downstream with it. Cavity 
surface currents and charges accompany the far field surface 


wave . 



V INTRODUCTION 

In a recent paper Schubert et al. [1973] took the first 
step in the development of an analytic theory for the asymmetric 
electromagnetic induction in the space defined by a spherical 
Moon and its downstream cylindrical cavity formed by the solar 
wind. Previous theories of lunar induction were spherically 
symmetric approximations [Schubert and Schwartz, 1969; Blank 
and Sill, 1969; Schwartz and Schubert, 1969; Schubert and 
Schwartz, 1972]. The asymmetric theory was developed in the 
quasistatic limit corresponding to large wavelength and low 
frequency of the oscillating interplanetary field. Thus Schubert 
et al . [1973] determined the asymmetric electromagnetic 

response of a two layer Moon model with an infinitely conducting 
core of radius b and an insulating layer of thickness (a-b) down- 
stream of which was a nonconducting cylindrical cavity. The 
frequency dependence of the response was treated qualitatively 
by the rough correspondence between frequency and the ratio of 
the core radius to the lunar radius (b/a). Small values of b/a 
correspond to low frequencies while larger values- of b/a correspond 
to higher frequencies. Despite the simplifications of this first 
asymmetric theory Smith et al. [1973] have obtained good agree- 
ment between theory and the experimental data of the Apollo 12 
Lunar Surface Magnetometer. 

In the quasistatic limit there are only two fundamental 
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orientations of the interplanetary magnetic field, parallel- ' 
and perpendicular to the lunar cavity axis. In the more complete 
time dependent scattering theory all directions for the plane 
wave propagation vector of the incident interplanetary field 
relative to the cavity axis must be considered. In the case of 
incident waves parallel to the cavity axis all linear polariza- 
tions are theoretically and physically equivalent since the in- 
cident electric and magnetic fields are both perpendicular to the 
cavity axis. However, when the propagation vector is not parallel 
to the cavity axis there are two distinct physical cases: the 
magnetic field may be perpendicular and the electric field may 
have a component parallel to the cavity axis, or the electric 
field may be perpendicular and the magnetic field may have a 
component parallel to the cavity axis. Thus there are three 
theoretical problems associated with the general solution for 
the electromagnetic scattering of a plane, linearly polarized 
wave incident on the asymmetric space defined by a spherical Moon 
and its downstream cavity. 

In this paper we obtain the general analytic solution for the 

propagation vector parallel to the cavity axis. 'We develop the 

theory in terms of an arbitrary Moon model subject only to the 

condition that the lunar electrical conductivity a , permittivity 

e and permeability p are functions of radius. The theory is also 

6 

limited to solar wind velocities less than 0(10 m/s) and to 
frequencies less than or equal to 10 Hz. The former condition is 
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of little consequence in^the lunar environment while the latter 
one encompasses the total frequency range of the Apollo surface 
and subsatellite magnetometer experiments. In future papers we 
will extend the theory to the case of the incident wave propa- 
gating perpendicular to the cavity axis. The asymmetric scatter- 
ing theory can then be used to invert Apollo surface and sub- 
satellite magnetometer data to yield lunar electrical conductivity 
distributions with the advantage that all magnetometer data, 
whether characteristic of the sunlit or night side lunar en- 
vironments, will be treated as a unified body of data to be 
inverted with a single theory. 
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^THEORETICAL MODEL 


In this section we lay the groundwork for the general 
analytic solution to the problem of electromagnetic scattering 
of a plane wave incident on a spherical Moon of arbitrary con- 

f 

ductivity o (r) and its downstream nonconducting cylindrical 
cavity when the incident wave propagation vector is parallel 
to the cavity axis. Following the notation of Schubert and 
Schwartz [1969} we represent the incident field by 
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where E and H are the electric and magnetic fields, a , a and 

— — • —X V:jr . 

a_ z are the unit vectors of a Cartesian coordinate system whose 
origin is at the center of the Moon, \ is the wavelength, cu 
is the circular frequency as measured by an observer fixed with 
the Moon, Hq is the amplitude of the magnetic field oscillation, 
v = Xu)/2tt and Ji is the magnetic permeability of free space. 
Throughout this paper MKS units are used. Also the time dependence 
e will henceforth often be understood. The wave propagates 

in the positive z-direction which is also the axis of the cylindri- 
cal downstream diamagnetic cavity. The Moon-cavity geometry and 
the geometry of the incident electromagnetic field are shown in 
Figure 1. Both spherical coordinates (r, 0, cp) and cylindrical 
coordinates (p, cp, z) will be used. The space outside the Moon 


and the cylindrical cavity contains only the incident inter- 
planetary electromagnetic field given by (1) . The fields 
associated with the current systems induced within the Moon, 
on the lunar sunlit hemisphere and on the cavity boundary cannot 
penetrate into this external region. The confinement condition 
on the front side is approximately correct since magnetohydrodynami 
waves cannot travel upstream in the supermagnetosonic solar wind 
plasma. Confinement to the interior of the cavity occurs because 
downstream travelling scattered waves in the supermagne to sonic 
solar wind are limited to the interior of the Mach cone extending 
downstream from the lunar limb. Since the Mach angle is small, 
the cylindrical surface is a good approximation to the Mach cone. 

The confinement is accomplished by currents which flow in 
the interplanetary plasma at the sunlit lunar surface and the 
cylindrical surface of the cavity. In the theoretical model 
these currents are assumed to be ideal surface current distri- 
butions, i.e. the current layers have negligible thickness. The 
boundary conditions are the continuity of tangential electric 
fields and normal magnetic fields on the sunlit lunar hemisphere 
and cavity boundary. 



THE MAGNETIQ FIELD IN THE LUNAR INTERIOR 

. t 

The electromagnetic field inside the Moon is most easily 
written in spherical coordinates. In general it consists of 
both transverse electric TE(E r =0) and transverse magnetic 
TM(H r = 0) modes, however, we will neglect the TM field. At 
the frequencies of interest f = w/2n < 1 Hz , the TM induction 
is negligible because of the relatively high resistivity 

near the lunar surface. 

On the sunlit lunar hemisphere, the radial component of 
the lunar magnetic field H r must be equal to that of the 
interplanetary field (the permeability of the Moon will be 
assumed to be that of free space) . The application of this 
boundary condition is facilitated by writing the TE part of the 
incident interplanetary field in multipole form.. Also, this 
spherical harmonic expansion of the TE incident field serves to 
indicate the form of the expansion for the lunar TE field. For 
the TE part of the incident magnetic field (1) we have [Schubert 
and Schwartz, 1972] 
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and (x) and P^(cos8) are the spherical Bessel functions and 
associated Legendre functions, respectively. 

The form of the TE magnetic field in the lunar interior is 
identical to that of incident field (2) except for the radially 
dependent functions. Following Schubert and Schwartz [1972] we 
write 
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and are coefficients to be determined. The (r) are 

solutions of the second order ordinary differential equation 
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subject to the conditions that are finite at the origin 

and G^(a) = 1 • The latter condition on is only a con- 

venient normalization condition; the coefficients A^ of each 
multipole contribution to H inside the Moon are as yet unknown. 
The propagation constant k is 

k = (co 2 pe + iuu no)'* ’ , 

where, for simplicity, M- and e have the free space values and 
o is an arbitrary function of r (it is not any more difficult 
to allow p and e to also be functions of r) . As an example, for 
the two layer Moon model with a perfectly conducting core of 
radius b and an insulating shell of thickness (a-b) 
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for frequencies of interest. 



THE CAVITY AS A CYLINDRICAL WAVEGUIDE BELOW CUTOFF 

Before we construct the general solution for the electro- 
magnetic field in the nonconducting cylindrical cavity downstream 
of the Moon, it is instructive to consider the physics of wave 
propagation in this cavity. Electromagnetic wave propagation 
in a hollow cylindrical waveguide is a classical problem dis- 
cussed, for example, in Stratton [1941]. A priori possible 

modes of propagation are transverse magnetic TM(H = 0) and transvers 

z 

electric TE(E z = 0). The TE and TM nomenclature in the cavity 
should not be confused with similar designations inside the Moon. 

The form of the incident field can be used to' restrict our 
consideration of the number of possible modes. In cylindrical 
coordinates the incident field (1) is 


{«} 


. , 2TTz -0Utx 

1 \ •i I 


H Q e 


incident 


A a 

flivcoscp p_ - pvsincp cck 
I sincp P_ + coscp cd J 


( 8 ) 


Thus we need only consider possible TE and TM modes of order 1, 
for which the lowest cutoff frequencies are 50.52 Hz and 105.15 Hz, 
respectively. The cutoff frequencies of the first order TE and 
TM modes are given by c/2na (c is the speed of light in vacuum) 
times the zeros of dJ^(x)/dx and (x) , respectively ( (x) 

is the Bessel function of first order) . 


Since the frequencies of interest are more than an order 

t 

of magnitude smaller than the lower of the above two cutoff 
frequencies, disturbances caused by the Moon cannot propagate 
down the cavity. Magnetic field perturbations in the cavity, 
due to the presence of the Moon, are oscillatory, with amplitudes 
which decay exponentially with distance down the cavity. . The 
first order TE disturbance decays essentially as exp(-1.84z/a) 
while the TM disturbance decays essentially according to exp (-3.8 
z/ a) . The electromagnetic field far downstream z-*” in the cavity 
is independent of the presence of the Moon . 
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THE ELECTROMAGNETIC FIELD IN THE 
FAR DOWNSTREAM CAVITY 


The electromagnetic field in the cylindrical cavity far 
downstream from the Moon can be obtained by a linear super- 
position of elementary cylindrical wave functions [see Stratton, 
'1941» pages 360-361]. From the form of the incident field (8) 
only first order TE and TM modes could contribute to the cavity 
field. Since there is no z-component in the incident electric 
field, continuity of the tangential electric field on the 
cylindrical boundary requires that the cavity field be a pure 
TE mode. Continuity of the radial magnetic field Hp on p=a 
uniquely determines the TE solution in the far downstream cavity. 
The condition that the tangential electric field component E^ 
be continuous on P=a is automatically satisfied. The solution 
is 
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E z = 0 » (10) 

where I^(x) is the modified Bessel function of order 1 and 
I ' (x) = dl-^(x)/dx. The field in the far downstream cavity (9) 
and (10) reduces to the incident field (8) in the limit ur» 0 or 
v-*c . . 

The TE cavity field (9) and (10) is the result of a surface 

wave on the cavity boundary induced by the incident field. This 

is most easily understood by considering the surface wave on a 

plane y=0 induced by a travelling wave (1) . In the region y>0, 

the field is everywhere the incident field and the half-space 

y<0 is vacuum. The field in the vacuum is independent of x , 

-UJt) 

proportional to e and dependent on y. The vacuum 

electric field has only an x-component which is given by 
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It is a surface wave travelling in the z-direction with speed v. 
Its amplitude is exponentially damped with distance from the 
surface into the vacuum. The attenuation depth is on the order 
of a wavelength. Note that the electric field is continuous on 
y=0. The vacuum magnetic field has a y-component given by 



which is continuous with the y-component of (1) on y=0. Since 
the normal magnetic field in the vacuum is attenuated with 
distance from the surface, v*H = 0 requires that the surface 
wave possess a z-component of magnetic field 



The surface wave is thus not a transverse wave, however, it is TE 
Clearly (9) and (10) represent waves on the cylindrical surface 
of the cavity which are induced by the passage of the inter- 
planetary field. The magnetic field components are 90° out of 

\ 

phase. 


The surface current density K on the cavity boundary o=a is 
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Note that K and K approach zero as u)-*0 or v-*c. The surface 
co z ^ 

charge density E on p=a is 



As in the case of the currents, E-*0 as iu-’O or v-*c. 

Figure 2 shows H /H~ and -E /uvH n as functions of p/a 

p u cp u 

with ^ - v 2 /c 2 as a parameter. The dependence of H q and 

on sincp. exp[i (— — ~ ujt) ] has been suppressed. The field components 

H and E have been normalized with respect to the amplitudes 
p cp 

of the corresponding interplanetary field components. Thus the 
continuity of H p and E^ on p=a makes the normalized field com- 
ponents unity on p=a, independent of the value of the wavelength 
parameter. When the wavelength parameter is zero, the normalized 
field components are identical to those of the incident field. 

As the wavelength parameter increases the field components of the 
surface wave are increasingly attenuated with distance from the 



cavity boundary. 




The normalized cavity field components H and E„ are shown 

co P 

2rra i 2 2 

in Figure 3 as functions of p/a with — — Jl-v z /c as a parameter. 

The cp, z and t dependences of H and E p have been suppressed. 

These field components are seen not to be continuous on o=a, 

except when the wavelength parameter approaches zero and in this 

limit H and E become uniform and equal to the incident field 
cp P . 

components. As the wavelength parameter increases, the amplitudes 
of these surface wave components at p=a decrease and the attenua- 
tion of these components with distance from the cavity boundary 
increases. At o=0, E p /uvH Q = H^/H q = - E^/uvHq = H p /H Q = O'. 

II 

Figure 4 shows the cavity field component z — 

i Hjl-vVc 2 

as a function of p/a with the same parameter as in Figures 2 
and 3. Again the cp, z and t dependences have been suppressed. 

Note that the incident field has no z-component. Also the H z 
in the cavity is 90° out of phase with all the other field 
components. For the wavelength parameter equal to zero, H 

z 

vanishes. As the parameter increases the normalized surface 
amplitude of H z increases, reaching a maximum value of about 1.1 
at a parameter value about 2.8, and then asymptotically decreases 
to unity. As with the other components, the attenuation of H z 
with distance from the cavity boundary increases with the wave- 
length parameter. Note that H z vanishes as p-*0. 
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The discontinuities in H and H at p=a are associated 

2 Cp 

with the cavity surface current densities K and K , respectively 

V z K 

The normalized surface current densities K /H^ and — 

z 0 


i H 


./l-v 2 /< 


are shown in Figure 5 as functions of the wavelength parameter 

(the dependences on z,t and cp are suppressed) . . A surface charge 

density £ is associated with the discontinuity in E^ at p=a. 

2 

The normalized surface charge density 2 c /vH Q is identical to 

K z /H 0 . As the wavelength parameter approaches zero the surface 

currents and charge approach zero. The normalized K and £ 

z 

monotonically approach unity as the parameter increases, while 
the normalized exhibits a maximum before eventually approach- 
ing unity. The magnitude of the azimuthal cavity surface currents 
exceeds that of the axial surface currents by a factor which 
increases monotonically as the parameter approaches zero. As in 

the case of H , K is 90° out of phase with K and £ . 
z cp z 


For a wavelength parameter smaller than 1/2 , the cavity 
field components H , H , E , and E are essentially the same 

D CO P Cp 

as the corresponding components of the interplanetary field. The 

parameter value 1/2 corresponds to a frequency of 0.02 Hz , for 

an assumed wave propagation speed of 400 km/sec. The wavelength 

parameter must be smaller by a factor of 2 (frequency smaller 

than 0.01 Hz) for H in the cavity to be negligible. 

z 



For hydromagnetic waves and discontinuities in the solar 
6 2 2 "5 

wind v < 10° m/sec and v /c < 10 . Using the smallness of 

2 2 

v /c , the expression (9) for the cavity magnetic field as 
z-*«> can be rewritten as 
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(17) 


In the limit of large wavelength i.e. 2Tra/\ « 1 , (17) is nearly 

a plane travelling wave with the magnetic field in the y-direction. 
With the aid of an identity from Cooke I" 1956] , (17) can be re- 

written in a form appropriate for satisfying matching conditions 
on the night side lunar hemisphere 
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MATCHING THE LUNAR AND FAR 
CAVITY MAGNETIC FIELDS 


The far cavity field satisfies the boundary conditions on 
the cylindrical cavity surface exactly. To match the lunar TE 
magnetic field to the far cavity magnetic field it is only 
necessary to add to the far cavity magnetic field the decaying 
TE mode of the cylindrical wavequide. The decaying TE mode 
satisfies = E^ = E z = 0 on p=a and consequently the sum 

of the far cavity field and the TE waveguide mode match the 
incident field correctly on the cavity surface P=a. The magnetic 
field of the decaying TE waveguide mode can be written [Stratton, 
1941, p. 541] 



1 

^ J 1 < V m £ a ,sinE 1 


? “hz 


1 ‘I 

• = H Y C e m * 

0 l-j in 

( J 1 (y m p / a )/ p ) c °sq) 


m=l 


Ih 

V. Z 


z 

( h J 1 ^ sinc P 

v m m 


where y for m=l,2,**** are the ascending zeros of J,' (x) , 
■ in J. 



im 




9 


(19) 


a 


( 20 ) 



and hi is the mth cutoff frequency cy m /a • As noted in our 
discussion of the cavity as a waveguide below cutoff, the 
frequencies of interest lie one or more orders of magnitude below 
the smallest cutoff frequency of about 50 Hz. Thus it is an 
excellent approximation to take h^ « y m /a an< ^ rewrite the 
magnetic field (19) as 

• f ” z/a . 

H=n H 0 l sin, C m e m • (2D 

m=l ' ' • 


This is identical to the quasistatic cavity disturbance field 
determined by Schubert et al. F1973]. To match with the lunar 


TE magnetic field on the night side hemisphere, (21) must be 
rewritten in spherical coordinates. Using Cooke [1956] we find 


H = v{-sinco H n 7 C Y 
— v. ,0 Z_i m Zj 

- m=l 1=1 


(- 1 + 1 ) ! 


P (cos A] 


Because we have used only the decaying TE waveguide mode 
to appropriately match the magnetic fields on the night side 
lunar surface, some mismatch of the electric fields is possible. 
This could be corrected by the inclusion of the decaying TM 
waveguide mode (and the lunar TM mode) “which would contribute 
a magnetic field a factor of ui /u> Smaller than the TE magnetic 
field. Thus for the frequencies of interest we make a negligible 
error by matching with only the decaying TE waveguide mode. 



COMPLETE ANALYTIC SOLUTION 

In the previous sections we developed representations 

for the TE magnetic field in the lunar interior (4) and in the 

cavity, the sum of (17) and (21), for an incident linearly 

polarized plane wave propagating in the positive z direction with 

velocity v«c (=3x10 m/s) , frequency f« 50 Hz and with the 

magnetic field vector in the y direction. These forms contain 

two infinite sets of constants A and C which must be evaluated 

n n 

by using appropriate boundary conditions. By the manner in which 

we constructed the cavity field, the boundary conditions at the 

cylindrical cavity surface, p=a , are automatically satisfied. 

The continuity of the normal component of the magnetic field at 

r=a over the entire Moon and the continuity of the tangential 

component of the magnetic field on r=a in the cavity, (K 8^ ^ , 

determine the doubly infinite set of constants. The algebraic 

detail of this procedure is given in the appendix, together with 

final formulae for the A and C . 

n n 
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SUMMARY 


In this paper we have, obtained the general analytic 
solution for the interaction of the Moon and its downstream 
cavity with a linearly polarized plane electromagnetic wave 
propagating in the direction parallel to the axis of the cavity. 

The solution is formulated in terms of a spherical Moon model 
with arbitrary radially dependent electromagnetic parameters 
and a nonconducting cylindrical downstream cavity. A number 
of approximations consistent with the physical circumstances 
of the Moon-solar wind interaction and our present ability to 
measure this interaction have been made. Outside the Moon and 
the cavity the electromagnetic field is the incident interplanetary 
field. The electrical conductivity of the near lunar surface 
is sufficiently low to suppress the TM lunar field. The incident 
wave velocity v is much less than the velocity of light in 
vacuum and the highest frequency of interest is of order 10 Hz. 

All but the first of these assumptions could be relaxed without 
presenting any special difficulty in carrying through the solution 
with the methods of this paper. * 

\ 

Of particular significance is the finding that for frequencie 
below about 50 Hz, the diamagnetic cavity is a cylindrical wave- 
guide below cutoff, i.e. cavity electromagnetic field perturbations 
due to the Moon decay exponentially with downstream distance. 


Thus the quasistatic representation of the cavity perturbation 
magnetic field due to the Moon given by Schubert et al. T1973] 
is valid for frequencies as high as 10 Hz. This is not to imply 
that the quasistatic solution for the lunar and cavity magnetic 
fields is valid at such high frequencies. Quite to the contrary, 
high frequency modifications in the lunar and cavity magnetic 
fields will be present at frequencies above about 0.01 Hz. 

Far downstream from the Moon, the interplanetary field 
forces a surface wave to propagate down the cavity boundary. 

The far cavity field, which is determined by this, propagating 
surface wave, is different from the interplanetary field, e.g. 
it has an axial magnetic field. component. These differences 
become more significant as the frequency of the field increases. 
The far cavity field is independent of the presence of the Moon 
whose cavity field perturbations will be damped within several 
lunar radii downstream. Differences between the interplanetary 
and cavity fields measured farther than a few hmmar radii down- 
stream must be attributed to the cavity itself,; they can yield 
no information about the lunar electrical conductivity. 
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APPENDIX 


To determine the unknown coefficients A and C first 

n n 

consider the continuity of the normal component of the magnetic 
field on r=a. We can write this condition in the form 
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On the lunar night side, within the cavity >fia0) and 
on r=a , both and H are continuous. It is sufficient to use 
only the continuity of H . This leads to the equation 
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To generate a set of equations for the unknown coefficients 
from (A-8) we note that on 0£9 ^tt/ 2 the odd associated Legendre 
functions form a complete set. Expanding (A-8) in terms of this 
set gives the infinite set of equations. 
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Equations (A-5) and (A-10) are sufficient for the determina- 
tion of the A^ and C f . After considerable algebraic manipulation 
a single set of equations can be derived for the , 
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for. l'=l, 3, 5, • • • In the limit v-*c and (A-ll) reduces 

exactly to (A-31) of Schubert et al . [1973] with the identifica- 
tion of u. as 1(1+1) a. . Equation (A-ll) represents an infinite 

v i/ 

sdt of linear equations (one for each l') for the infinite set of 
unknowns C m . Since the C m eventually decrease in importance as 
m increases the set of equations can be solved numerically by 
truncating in m and t ' . This was done in the quasistatic limit 

o * 

by Schubert et al . ("19731 and will be done for various Moon models 
using the present theory in a future paper. Once the have been 
found, the A^ may be obtained from (A-5) . After some algebraic 
manipulation we find 
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Again, in the limit of v-*c and X"* 00 , the equations for the 
quantities P^A^ formally reduce to equations (A. SO) and (A. 31) of 

t 

Schubert et al . [1973] . 


FIGURE CAPTIONS 


Figure 1. The geometry of the asymmetric lunar electromagnetic 
scattering problem. 

Figure 2. The far downstream cavity field components H q and 
normalized to the corresponding incident field 
components as a function of normalized distance from 
the cavity axis with — as a parameter. 

Figure 3. The far downstream cavity field components H and E q 

normalized to the corresponding incident field components 
as a function of normalized distance from the cavity 
axis with - ^ a ^l-v^/c^ as a parameter. 

Figure 4. The far downstream cavity axial magnetic field 

j 2 2 

H /i H^/l-v, /c as a function of normalized distance 
z u* 

from the cavity axis with as a parameter. 

The dependence of H on z, cp and t has been suppressed. 

z 

^ i 

\ 

Normalized surface current density and charge 
density on the far downstream cavity boundary as a 
function of ,/l-v^/c^ . The dependence on z, co 

and t has been suppressed 


Figure 5. 
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